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ABSTRACT 

We measure the intracluster medium temperature distributions for 62 galaxy clusters in the HI- 
FLUGCS, an X-ray flux-limited sample, with available X-ray data from XMM-Newton. We search for 
correlations between the width of the temperature distributions and other cluster properties, including 
median cluster temperature, luminosity, size, presence of a cool core, AGN activity, and dynamical 
state. We use a Markov Chain Monte Carlo analysis which models the ICM as a collection of X-ray 
emitting smoothed particles of plasma. Each smoothed particle is given its own set of parameters, in- 
cluding temperature, spatial position, redshift, size, and emission measure. This allows us to measure 
the width of the temperature distribution, median temperature, and total emission measure of each 
cluster. We find that none of the clusters have a temperature width consistent with isothermality. 
Counterintuitively, we also find that the temperature distribution widths of disturbed, non-cool-core, 
and AGN-free clusters tend to be wider than in other clusters. A linear fit to (t^t — kTmed finds 
CfeT ~ 0-20kTmed + 1.08, with an estimated intrinsic scatter of ^ 0.55 keV, demonstrating a large 
range in ICM thermal histories. 

Subject headings: Galaxies: clusters: general — Galaxies: clusters: intracluster medium — X-rays: 
galaxies: clusters 



1. INTRODUCTION 

Galaxy clusters have proven to be invaluable tools for 
cosmological studies. The largest gravitationally bound 
objects in the Universe, clusters trace their origins to 
fluctuations in the primordial density field. As a re- 
sult, the evolution of galaxy clusters is intimately linked 
with the growth of large scale structure. Specifically, the 
cluster mass function, n{M,z), is sensitive to flM and 
f^A, the matter and dark energy densities of the uni- 
verse, and as, the mass fluctuation amplitude on scales 
of ~ 8h~^ Mpc. Thus measuring the mass of clusters 
allows us to probe large scale structure in the universe. 
Mass measurements for a range of redshifts can then help 
to trace the formation and evolution of this structure. 
The viability of using clusters as cosmological probes 
relies on obtaining accurate cluster masses. Typically, 
this requires using an observable mass proxy, such as X- 
ray temperature or luminosity, which scales predictably 
with total cluster mass. However, significant scatter in 
such scaling relations is a major source of error in clus- 
ter mass measurements. It has been found that cluster 
X-ray morphology and temperature substructure, par- 
ticularly the presence of dense, cool cores, signi ficantly 
affect s the scatter in X-ray scali n g relations (.Markevitchl 
TOOa lArnaud fc EvrardI [19991: lAlien fc FabianI iToM 
Zhang et all 120071 : IVentimidia et al.l I2008D . A better 
grasp of cluster physics is clearly essential to understand- 
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ing the scaling relations and thus obtaining more precise 
cosmological measurements. 

The complex physics of galaxy clusters results in a 
wealth of temperature substructure. Cluster-cluster in- 
teractions, including mergers, radiative cooling, and 
AGN f eedbac k can induce shocks and cold fronts 
(|Markev itch fc Vikhlininll2007[ ). dense c ool cores. X-ray 
cavit i es, and soft X-ray filaments (e.g., iCrawford et al.l 
[200aiFabian et alllMl . Additionally, the cold gas ex- 
pected from traditional cooling flow models is largely 
missing from cluster cores, observations instead reveal- 
ing a cut off in the teinperat ure distributions at T '--^ 
Tmax/'i (Peter son et al.l [20031 ) . despite some evidence of 
minor cooling flow ac tivity in the form of star-formation 
(iBildfell et al.l 120081) . Ha filaments (jMcDonald et all 
12010! ) , and a small amoun t of cold X-ray gas (e.g ., 
Sand ers et aXE OQ^. 200 9lbl . IMol : IPeterson et al.l[2003l) . 
Thus it appears that while cooling flows are not com- 
pletely quenched, some amount of heating is neces- 
sary to explain the missing low-temperature gas in clus- 
ter cores. There are many proposed heating mecha- 
nisms, some of the most common being turbulence, ther- 
mal conduction, or some form of AGN feedback. The 
most popular is AGN feedb ack via inflation of cavities 
(jMcNamiara fc NulsenI |2007[ ). It has been shown that 
there is enough energy available in X-ray cavities to ad- 
equately balance r adiative cooling ([Raffertv et al.l 120061 : 
IBirzan et al.ll2004 ). but the exact nature of how this en- 
ergy is distributed throughout the ICM is unknown. 
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It is clear that cluster temperature structure is strongly 
tied to the dynamical history of the cluster. Thus, a 
more thorough understanding of cluster temperature dis- 
tributions may lead to reduced scatter in the X-ray scal- 
ing relations used for cosmology measurements and shed 
light on the cooling flow problem, cluster heating, and 
cluster-cluster interactions. Most measurements of clus- 
ter temperature distributions treat the ICM as a mul- 
tiphase plasma, with anywhere from one (isothermal) 
to four temperature components. Typically, this also 
involves assumptions on the spatial distribution of the 
multi-temperature plasma. Using the S moothed Parti- 
cle Inf erence (SPI) method introduced bv lPeterson et al.l 
(|2007D it is possible to avoid such spectral and spa- 
tial assumptions by treating the ICM as a collection of 
plasma parcels, each emitting X-rays according to an 
independent set of parameters, including emission mea- 
sure, temperature, spatial position, and size. These pho- 
tons are then propagated through the XMM-Newton in- 
strument response and the model parameters adjusted 
via Markov Chain Monte Carlo. The resulting distri- 
bution of smoothed particle parameters provides a good 
description of the cluster. In this way, we are able to mea- 
sure the full temperature distribution, across all temper- 
atures, while at the same time allowing for both temper- 
ature and luminosity substructure and asymmetry. We 
aim to characterize cluster temperature distributions in 
several ways. First, we investigate how well we are able to 
recover the shape of the temperature distribution. Sec- 
ond, we measure the width of the distribution, <TkT, as 
well as the cluster median temperature, total emission 
measure, and r25oo- Third, we search for correlations of 
(TfeT with other cluster properties, such as median tem- 
perature. Information on each cluster in our sample is 
gathered from the literature in order to determine if clus- 
ter dynamical state, central AGN activity, or cool-core 
status has any effect on the temperature distributions. 

Assumed cosmology throughout the paper is Hq = 70 
km s-iMpc-\ flM = 0.3, and = 0.7. 

2. DATA 
2.1. Data Sample 

Our cluster sample is the Highest X-ray FLUx 
Galaxy Cluster Sample ( HIFLUGCS) compiled by 
IReiprich fc Bohringed (|2002[ ). HIFLUGCS is a statisti- 
cally complete flux-limited sample, containing the 63 X- 
ray brightest clusters in the sky (excluding the Galactic 
plane). The HIFLUGCS also contains a variety of clus- 
ter morphologies, including merging clusters, cool core 
and non-cool core clusters, and both galaxy clusters and 
galaxy groups. The clusters have redshifts in the range 
0.003 ^ z < 0.2 and span a wide range in temperature, 
0.8 keV < Tx < 13 keV. We selected this sample be- 
cause we wanted to both span the range of cluster types 
and have enough photons in the data to allow detailed 
measurements . 

Of the 63 clusters in the sample, 62 have been 
observed with XMM-Newton (ah except Abell 2244). 
For clusters with multiple available observations, we 
choose the longest observation (after flare screening) 
that is nearest the cluster center. We use the clus- 
ter optical redshifts as given in the NASA/IPAC Ex- 
tragalactic Database (NED). Galactic hydrogen col- 



umn densities are inferred from t he 21c m radio mea- 
surements of [Pickev &: LockmanI (|1990f) . as given in 
IReiprich fc Bohringeii ( 20021 ). The observation details 
are given in Table [T] 

2.2. Data Reduction 

Our analysis requires only a photon event list and 
exposure map for each XMM-Newton detector (MOSl, 
M0S2, pn, RGSl, and RGS2). These were created us- 
ing the SAS 10.0 pipeline tasks emchain, epchain, and 
rgsproc. Light curves were created from the resulting 
event files to identify periods affected by soft proton 
flares, which were then removed from the data for the 
remainder of the analysis. Net exposure times, after 
flare screening, are listed in Table [1] for each observa- 
tion. In addition, EPIC event files were filtered to include 
only photon event patterns — 12 (to exclude non- X-ray 
events) and photon energies in t he range 0.3 — 10.0 keV 
(MOS) and 1.1 - 10.0 keV (pn) ()Andersson &: Madeiskil 
I2OOI . 

3. ANALYSIS 

We employ the smoothed particle inference method 
presented in Peterson, Marshall, and Andersson (2007) 
as an alternative to traditional X-ray analysis procedures 
for galaxy clusters. The smoothed particle method mod- 
els the intracluster medium as a collection of X-ray emit- 
ting parcels of plasma. We refer to these parcels as 
smoothed particles, and each represents a large volume of 
plasma in the ICM. Each smoothed particle is character- 
ized by a set of parameters, including temperature, emis- 
sion measure, gaussian spatial width, redshift, and abun- 
dance. The smoothed particles emit photons according 
to the chosen spectral and spatial model f t|3.ip . which 
are then propagated through the XMM-Newton instru- 
ment response. The resulting simulated data are then 
iteratively compared with the XMM-Newton data via 
a Markov Chain Monte Carlo (MCMC) process {MM- 
The EPIC and RGS data are processed separately for 
each observation and the results combined prior to mea- 
suring cluster properties (t )3.3.ip . We use the combined 
MCMC results to measure the cluster radius, emission 
measure, median temperature, and temperature distri- 
bution width, (7kT (J 



3.1. Smoothed Particle Model 

The model for each smoothed particle consists of mul- 
tiple components: the ICM model, a Galactic and extra- 
galactic X-ray background, and the instrumental back- 
ground. Though some parameters are global (the same 
for all smoothed particles), the four-component model 
applies independently to each smoothed particle. The 
fractions of photons going to each model component, as 
well as the relative normalizations, are free parameters. 
The model parameters are summarized in Table [H 

3.1.1. ICM Model 

For our analysis, the smoothed particles are spa- 
tially modeled as gaussians. The gaussian width as 
well as the two-dimensional position of each smoothed 
particle is allowed to vary independently. The ICM 
spectrum of each smoothed particle is p roduced by 
a WABS ((Morrison fc McCammon. I1983D absorbed. 
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isoth e rmal MEKAL (iMewe et all 119851 119861: iKaastral 
119921 iLiedahl et aIlll995D model Individual smoothed 
particles thus have a single temperature which may be 
different from every other smoothed particle. The prior 
for the distribution of smoothed particle temperatures 
is logarithmically uniform between 0.1 keV < kT < 
10.0 keV, except for EPIC observations of clusters with 
kTavg > 4 keV. Since these clusters may have some gas at 
temperatures > 10 keV, which EPIC can detect, the prior 
for these clusters is set to 0.1 keV < kT < 19.5 keV. The 
absorbing column density, nn, and the smoothed particle 
metallicity, Z , are also allowed to vary; however they are 
global parameters. Galactic column densities and red- 
shifts are those listed in Table [TJ The smoothed particle 
metallicities are with resp ect to the solar abundances of 
lAnders fc Grevessd (|1989D and use a uniform prior from 
O.OZq -2.OZ0. 

3.1.2. X-ray Background 

The XMM-Newton data include not only photons 
emitted by the ICM, but also photons from the X-ray 
background. The X-ray background has three compo- 
nents: a soft Galactic X-ray background, a hard extra- 
galactic X-ray background, and the instrumental back- 
ground. To account for these background photons, addi- 
tional components are added to the model described in 
t jB.l.ll The normalization of each component is a free 
parameter. 

The Galactic X-ray background consists primarily of 
soft X-ray photons of energies < 1 keV with a spec- 
trum consisting of several therma l components from 
0.05 - 0.5 keV dLu mb et^ [200l IKuntz fc SnowdenI 
20001: ISnowden et a l. 2008). Following the model of 



Andersson et al.l (|2009l ). based on analysis of blank sky 



data files, we use a spatially uniform, unabsorbed, 
isothermal MEKAL model to describe the Galactic X- 
ray emission. For all smoothed particles, the soft X-ray 
background temperature is fixed at 0.15 keV, the red- 
shift is fixed at z = 0, and the metallicity is fixed at 
Z = I.OZq. 

Unresolved extragalactic sources, such as AGN, are re- 
sponsible for hard (energy > 1 keV) X-ray background 
emission. The spectrum is well described by a power-law 
with an index of P = 1.47 and is spatially uniform. A 
model component with a spatially uniform, WABS ab- 
sorbed, power law spectrum (with P = 1.47) is used to 
represent the extragalactic background. The absorbing 
column density is fixed to the value of nn for each cluster 
(Tabled. 

The R GS and EPIC instrume ntal backgrounds (de- 
scribed in lden Herder et al.|[200TI and lLumb et al.ll200a 

respectively), consist of several components. A soft pro- 
ton background is modeled as a power law with P = 
-0.205 (-0.35) for EPIC (RGS). Electronic noise is mod- 
eled as an exponential, F oc e~^/^', with E, = 158 (500) 
eV for EPIC (RGS). We also model instrumental line 
emission. In addition, RGS has four calibration sources, 
which are modeled as 2D spatial gaussians. For further 
discussion on our implementa tion of the RGS and EPIC 
instrumental backgroun ds, see lPeterson et al.l (2003) and 
lAndersson et al.l ()2007f ). The parameters for the instru- 
mental background component are fixed for all smoothed 
particles, and are the same for all observations. They 
only differ between the EPIC and RGS instruments. 



3.2. Smoothed Particle Inference 

The four spectral and spatial model components de- 
scribed in ij3.1l are combined for the MCMC analysis. 
For each smoothed particle, photons are emitted accord- 
ing to this model, which are then propagated through 
the instrument response (either EPIC or RGS) and com- 
pared with the input XMM-Newton data. 

Some photons in the data are removed based on spatial 
position prior to beginning the MCMC analysis. The in- 
put images are 30 x 30 square arcminutcs for EPIC and 
5 arcminutes by 30 arcminutes for RGS. EPIC photons 
outside a 20' x 20' square centered on the cluster are re- 
moved, or outside a 30' x 30' region for clusters in which 
the core radius, as given inlReiprich & Bohringcr (2002), 
extends farther than 10' from the observation center 
(only applies to 4 clusters). RGS photons are removed 
if more than 3' (2.5') from the center in the dispersion 
(cross-dispersion) direction. Photons are also removed 
according to energy. EPIC MOSl and M0S2 photons 
are only kept if they have energies 0.3 keV < E < 10.0 
keV. Similarly, pn photons in the energy range 1.1 keV < 
E < 10.0 keV are kept (cf. I Andersson fc Made iski 2004). 
The number of EPIC photons is increased by less than 
2% if the pn energy range is extended down to 0.3 keV, 
and thus the difference in MOS and pn energy cuts does 
not have a significant affect on analysis. Wavelength cuts 
for RGS are 5A < A < 65A (0.2 keV - 2.5 keV). There 
are several observations in which the cluster is outside of 
the RGS field of view and therefore only the EPIC data 
are used for these observations. 

The number of smoothed particles (Nsp) and the over- 
sim ratio (r, the ratio of simulated photons to data pho- 
tons) are chosen such that 



Npr 
iNsp 



lo^ 



(1) 



where Np is the number of photons in the data, after 
making the cuts discussed in the above paragraph. The 
reasoning for this criteria is described in ^'SA\ r is typi- 
cally chosen to be 10. However, to mitigate the already 
large computing time, if Np > 1.5 x 10^, we use r = 1 
and 10 times fewer smoothed particles, in accordance 
with Eqn. [TJ We also set a lower limit to the number of 
smoothed particles at Nsp > 10. This limit comes into 
effect for observations with very few photons (typically 
RGS observations). In these cases, r = 100 and Ngp is 
increased accordingly. Values of Ngp range from 10 to 
several hundred. 

The MCMC typically converges by iteration 500; how- 
ever to be conservative, we automatically disregard all 
iterations prior to 1000. Convergence is verified by bin- 
ning the iterations (300 iterations per bin) and measur- 
ing the average ^^.d primary cluster parameters (as 
described in i i3.3p . except cluster size, for each bin. The 
chain is considered to be converged when these values 
agree within Icr for 10 neighboring bins (the average 
must be within 2%). The MCMC is allowed to run until 
a minimum number of iterations (after convergence), N, 
is met. This minimum depends on N^p for the observa- 
tion, such that the total number of smoothed particles 
to be used in our analysis, NgpN, is 3 x 10^ for EPIC 
and 1 X 10^ for RGS. However, a minimum of 2000 itera- 
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tions are required for all observations, regardless of N^p. 
For high quality observations with many smoothed par- 
ticles, this results in using iterations 1000 to 3000 for our 
analysis. Computation times vary from several days to 
a few months. We run the MCMC on a pool of parallel 
machines. 

3.3. Measurement of Cluster Properties 

The desired cluster properties can be calculated di- 
rectly from the MCMC results by analyzing the distri- 
bution of individual smoothed particle parameters. The 
primary properties of interest are median cluster temper- 
ature {kTjned), the width of the temperature distribution 
(cfcr), a characteristic cluster radius (?'25oo), and the to- 
tal emission measure (EM) or luminosity. Uncertainties 
for most of these properties are also relatively simple to 
calculate using the standard deviations of the MCMC 
resuhs (g3Xni). 

3.3.1. Combining EPIC and RGS 

To take full advantage of the different energy sensitiv- 
ities of RGS and EPIC, the MCMC resuhs from EPIC 
and RGS observations are combined for the purpose of 
measuring cluster properties. In fact, to measure the en- 
tire X-ray temperature distribution, the use of both is 
essential. However, the EPIC field of view is larger than 
that of RGS (the RGS field of view is a 5' strip across 
the center of the EPIC field of view), and so care must 
be taken when combining the smoothed particle param- 
eter distributions. Within the RGS field of view, where 
there are smoothed particles from both the RGS and 
EPIC MCMC results, the aim is to give more weight to 
the RGS observation at lower temperatures and to the 
EPIC observation at higher temperatures. The RGS in- 
strument determines temperature primarily via K- and 
L-shell emission lines of C, N, and O, as well as Fe L- 
and M-shell emission lines. EPIC is primarily sensitive 
to Fe K-shell emission lines. The temperature sensitivity 
of RGS and EPIC can therefore be roughly characterized 
by the ionization fraction of Fe L and Fe K, respectively, 
as a function of temperature. Note that the shape of 
the bremsstrahlung continuum is important when using 
EPIC to measure a single cluster temperature. However, 
it is difficult to disentangle multiple temperature compo- 
nents based on the continuum; consequently, it is more 
appropriate to use the Fe K emission lines to characterize 
the ability of EPIC to measure the temperature distri- 
bution. In the overlapping region, we therefore weight 
the RGS results (the individual smoothed particle pa- 
rameters) by the Fe L ionization fractions and the EPIC 
results by the Fe K ionization fractions (see Figure [1]). 
This ensures that where RGS data is available, it is used 
preferentially for determining the temperature distribu- 
tion where it is more sensitive than EPIC, namely at 
temperatures < 1.76 keV, and similarly, that EPIC data 
is always preferred at higher temperatures. Outside of 
the RGS field of view, the EPIC results are weighted by 
the Fe K -|- Fe L ionization fractions. In this way we 
are able to take advantage of the greater sensitivity of 
RGS at low temperatures, while at the same time utiliz- 
ing the higher temperature sensitivity and larger field of 
view of EPIC. Note that RGS data contains very little 
spatial information, other than that based on its field of 



view, and therefore this method should not be applied if 
investigating the spatial dependence of the temperature 
distribution, e.g. for making temperature maps. For the 
purpose of this work, it is only relevant whether or not 
RGS detects gas at a given temperature, not where this 
gas is located. 

3.3.2. r-2500 

The characteristic radius of each cluster is chosen to 
be ^2500, the radius within which the mean mass density 
is 2500 times the critical density of the universe at that 
redshift. To dete rmine r25oo, we emp loy the M — Tx 
scaling relation of lArnaud et al.l ()2005D . 



h{z)M. 



2500 



1.69 



kT 
5 keV 



1.70 



(2) 



and calculate the radius within which this is true for the 
measured cluster temperature. 



'^2500 



3 * Af. 



2500 



1/3 



(3) 



The temperature used in lArnaud et al.l (|2005[ ) is the core- 
excised spectroscopic temperature from 0.1 — 0.5r2oo- 
However, for the purposes of this work, it is not crucial 
to have a precise estimate of r25oo, and we use here the 
emission-weighted median temperature of all smoothed 
particles (all of which are within the spatial cuts ii3.2p . 
Excising t he core prior to measu ring the temperature, as 
is done in lArnaud et al.l ()2005l ) , only changes our r25oo 
measurements by a few percent. r25oo is the first cluster 
property measured, and all cluster properties are mea- 
sured within this radius by only considering smoothed 
particle emission within r25oo • 

3.3.3. Median Temperature 

There are many different ways to define a single tem- 
perature for a cluster. We choose to use a temperature 
kTmed defined by the median log kT of the smoothed 
particles. Smoothed particles from both RGS and EPIC, 
within r2500 7 are used. For each smoothed particle, logkT 
is weighted by both its emission measure and Wrgs(^T) 
or W^picikT) (where kT is the smoothed particle's tem- 
perature and W^Gs{kT) or W^pidkT) is the ionization 
fraction weight for a smoothed particle of temperature 
kT). Recall that for EPIC smoothed particles which 
are outside the RGS field of view, Wepic is equal to the 
sum of both the Fe K and Fe L ionization fractions (Fig- 
ure [T]) , regardless of the smoothed particle temperature. 
The median of this weighted distribution then defines the 
cluster's median temperature, kTmed- 

3.3.4. Width of Temperature Distribution 

All of the other primary cluster properties - physical 
size, total emission measure, and temperature - can be 
found using traditional analysis methods. The main ad- 
vantage of our method is that by constructing a contin- 
uous temperature distribution we can also measure the 
width of this distribution, UkT- 

Both RGS and EPIC smoothed particles (again within 
''2500) are used in the determination of (Jut- As with the 
kTmed measurement, each smoothed particle's tempera- 
ture is weighted by its emission measure and W^Gs{kT) 
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or WEPic{kT). The standard deviation of the weighted 
temperature distribution is then taken to be Ufcr, as given 
in Eqns. |4]and[5] Here, kTj and EMj are the tempera- 
ture and emission measure of the j*'* smoothed particle. 
Wj is either WncsikTj) or WEmcikTj). The sum is over 
ah RGS and EPIC smoothed particles. This value is then 
corrected for the isothermal spread function, as described 



Nrgs+^epic 

[kTj-Wf EM J Wj 



2 



N 



(4) 



kT ■ 



j 

Nrgs+I^epic 

^ kT,EM,W, 

i 

N 

Y^EM.W, 



(5) 



3.3.5. Emission Measure 



Rather than explicitly calculating the luminosity, we 
choose to use a closely related quantity, emission mea- 
sure, given by Eqn. [51 



EM ^ / neUidV 



(6) 



It is directly related to the X-ray luminosity through the 
cooling function, A(r). For most clusters, A ~ T^/^, 
i.e. only bremsstrahlung emission is significant, but for 
lower temperature clusters the effect of line emission be- 
gins to alter the cooling function. Therefore when com- 
paring clusters over a wide range of temperatures, where 
the cooling function is not constant, it seems more ap- 
propriate to use the emission measure rather than the 
luminosity. 

It is relatively straightforward to calculate the emis- 
sion measure for each smoothed particle (and thus also 
for the entire cluster) from the MCMC results. One of 
the smoothed particle parameters is the MEKAL normal- 
ization f §3.1.1|) . This is related directly to the emission 
measure (Eqn. [T]). 



10- 



norniMEKAL = 



4TrD\{l + zy 



EM 



(7) 



The total emission measure for the cluster is simply the 
weighted sum of the smoothed particle emission mea- 
sures, averaged over the number of iterations used, as in 
Eqn. [3 



-'Vrgs 



EMi 



total 



- J2 WkEMk 

' k 

(8) 

where the sums are over all the RGS (j) and EPIC (fc) 
smoothed particles, and TVhgs and A^epic are the num- 
ber of RGS and EPIC iterations. The sum over RGS 
smoothed particles yields the contribution to the total 
emission measure of cool gas within the RGS field of 



view. The sum over EPIC smoothed particles comprises 
the contributions of the higher temperature gas within 
the RGS field of view and gas of all temperatures out- 
side the RGS field of view. 

3.3.6. Statistical Uncertainties 

One of the advantages of using a Markov Chain Monte 
Carlo is that it comes with a built-in estimate of statis- 
tical uncertainties. Once the chain has converged, the 
parameters from any individual iteration comprise a sta- 
tistically acceptable representation of the cluster. The 
(weighted) contribution of RGS to EM , kTmed, and akT 
is calculated separately for each RGS iteration and the 
standard deviation of these values then comprises the 
RGS contribution to the error. The same is done for 
EPIC, and error on the combined quantity is then cal- 
culated using standard propagation of error techniques. 
In the case of ^2500, the error is calculated using propa- 
gation of error on Equations [2] and [3l including the er- 
rors on kTmed and t he two M — Tx fit parameters from 
lArnaud et all (f2005l) . 

3.4. Isothermal Spread Function 

There is an intrinsic limit to how well the tempera- 
ture distribution can be resolved. This limit is affected 
by two major factors. The first is the statistical ability 
of the data to determine the temperature distribution, 
given the imperfect instrument response and finite num- 
ber of photons. Secondly, even with a perfect instrument 
and infinite number of photons, the spectroscopic infor- 
mation from atomic transitions has an incomplete set of 
information about the temperature distribution. The ef- 
fect of this limit is to introduce extra broadening to the 
temperature distribution, similar to the effect of a point- 
spread function, so that the measured temperature width 
of even an isothermal cluster will be nonzero. We must 
correct for this 'Isothermal Spread Function' (ISF) in 
order to obtain a measurement of the cluster's intrinsic 
temperature width. We quantify the ISF by simulating 
a set of isothermal clusters, whose true akT = 0. The 
simulated clusters are analyzed as described in §3.31 and 
their (JkT values are taken to be measurements of the 
ISF, aisF- 

Clusters were simulated according to an isothermal 
MEKAL spectral model. The spatial distribution of each 
cluster is described by a /3-model, with = 60" and 
/3 = 0.8. po is chosen to get the desired number of pho- 
tons, Np. As a measure of the ability to measure ICM 
temperatures, the ISF is a spectral phenomenon, and is 
therefore not expected to be sensitive to spatial charac- 
teristics of the X-ray emission (e.g. Tc or f3). Energy and 
spatial cuts are the same as in §3.21 We determined the 
effect of four parameters on the ISF: the cluster temper- 
ature (fcr), the number of smoothed particles (Nsp), the 
number of photons in the data (Np), and the ratio of sim- 
ulated photons to photons in the data (the 'over-simulate 
factor' r). Combining these effects, we find the ISF can 
be approximated as aisp ~ ai{Npr / N ap)a2{kT) . 

3.4.1. Np, Nsp, r 

The first group of simulated clusters determined the ef- 
fect of Np, Nsp, and r on the ISF. The temperature of all 
simulated clusters in this group was set to 1 keV. The val- 
ues of Np, Nsp, and r for the model were chosen from the 
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following: Np = {10^, 10"*, 10^}. Nsp = {10,100,1000}, 
and r = {1, 10, 100}. All possible combinations of these 
parameters were simulated for both RGS and EPIC (with 



the exception of the RGS r = 100, Np = 10^ TV^. 



1000 



and EPIC r = 100, Np = 10^, Nsp = 10 simulations, due 
to the long computation time required), resulting in 26 
measurements of the ISF for EPIC and 26 for RGS. 

We found that the ISF depends on the ratio Npr/Nsp, 
which is the average number of simulated photons per 
smoothed particle (Figure [5]) . This ratio ranged from 1 
to 10^ in our simulations. The best fit as a function of 
Npr/Nsp for RGS is 



ai(^)a2(lkeV) 

sp 

and for EPIC is 



Npr 

Nsr, 



Nr,r 



N r 

ai(-r^)a2(lkeV)=5.50 



-0.43 



-0.49 



0.05, 



(9) 



0.09. (10) 



Unsurprisingly, the more simulated photons per 
smoothed particle, the better the temperature distribu- 
tion can be constrained. There proved to be too few 
simulated photons per smoothed particle in the r — 1, 
Np = 10"^, and Ngp = 1000 simulations to constrain the 
temperature distribution (these were excluded from the 
aiSF fits for this reason). Thus, it appears we must have 
Npr/ Nsp > 1, and given the swift increase in the size of 
the ISF, 10^ < Npr/Nsp < lO"* is preferable. Higher val- 
ues may decrease the ISF somewhat; however, they also 
drastically increase the already large computation time. 
We therefore choose Nsp and r for each dataset such that 
Npr/Nsp ~ 10*. Thus we are able to minimize that con- 
tribution to the ISF of the number of photons in the 
data by a choosing an appropriate number of smoothed 
particles and over-simulate factor. 

3.4.2. kT 

We minimize the contribution to the ISF due to the 
number of photons in the as described above. The cluster 
temperature is thus primarily responsible for determining 
the ISF. We therefore simulated clusters with tempera- 
tures (in keV) of kT = {0.5, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 13}. 
For each temperature, we simulated the same cluster 
(with a given po, rc, and /3) three times, with the fol- 
lowing instruments and temperature priors: 1) RGS, 
0.1 keV <kT < 10.0 keV, 2) EPIC, 0.1 keV < kT < 10.0 
keV, and 3) EPIC 0.1 keV < kT < 19.5 keV. po was cho- 
sen at each temperature such that Np = 10^ for EPIC. 
Nsp and r were then chosen individually for each simula- 
tion according to i j3.21 such that Npr/Nsp — lO"*. Each 
EPIC simulation was then combined with the RGS sim- 
ulation of the same kT and cr/s f measured according to 
i j3.3.4l The resulting (t/sf as a function of cluster tem- 
perature, measured in keV, for the EPIC temperature 
prior of 0.1 keV < kT < 10.0 keV is 

<JlSF{kT) = 0.41fcT2-13g-0.35fcT^ 

and for the 0.1 keV < kT < 19.5 keV temperature prior 
is 

aisFikT) = 0A3kT^°°e-°-^^''^. (12) 
Both functions are shown in Figure |3l 



3.4.3. ISF Correction 

For the MCMC analysis, the Npr/Nsp contribution to 
the ISF if minimized by setting r and Nsp such that 
Npr/Nsp ~ 10*, where Np is estimated as one quarter 
of the total photons in the data (since there are four 
model components, and only one component represents 
the ICM). The ISF is then calculated using Eqn. [TTjor 
1121 For the temperature, we use the combined median 
temperature as described in ij3.3.3l The ISF is subtracted 
from the raw akT ( ^3.3.4p according to Eqn. [131 All 
reported akT are corrected for the ISF in this manner. 
An additional systematic error term is incorporated in 
the reported <7kT errors to account for uncertainties in 
the ISF measurements and any residual error due to the 
number of photons. 



kT,raw 



4. RESULTS 



'ISF 



(13) 



4.1. Cluster Morphology and Dynamical State 

In order to investigate possible connections between 
our results and cluster morphology or dynamical state, 
we have taken advantage of the well-studied nature of 
the HIFLUGCS sample and found information on each 
cluster's X-ray and radio morphology in the literature 
(see references in Table [3]). For each cluster, its ra- 
dio features are classi fied according to the taxonomy of 
IKempner et"all ()2004D . the major categories being radio 
features that are associated with AGN and features that 
are associated with the ICM, which are usually indica- 
tive of merging activity (e.g. giant radio halos). X-ray 
results for each cluster are used to determine if the cluster 
is known to have X-ray cavities (with or without coin- 
cident radio emission) and whether it is known to have 
a disturbed morphology (such as filamentary structures, 
secondary X-ray brightness p eaks, or significant iso pho- 
tal centroid shifts, e.g. as in IVikhlinin et all (12001) 1 or 
a regular, symmetric X-ray morphology. It is also noted 
if the cluster is known to be undergoing a merger based 
on detailed X-ray observations. 

Based on this information, we divide the cluster sam- 
ple two different ways. First, we have separated clus- 
ters which are known to be relaxed from those which 
are disturbed in any way. A cluster is designated as dis- 
turbed if it has any of the following features: evidence 
of merging activity (from detailed X-ray observations or 
large-scale radio emission), or a disturbed X-ray bright- 
ness morphology. A cluster is designated as relaxed only 
if it has none of the latter and has a regular, symmet- 
ric X-ray morphology (which may include cavities). A 
cluster typically satisfies several of the criteria for its 
designated category. For all but six cases, this classi- 
fication matches the disturbed / un disturbed designations 
of t he HIFLUGC S clust ers used in lVikhlinin eFall (l2009l) 
and iZhang et al.l (|2011[ ) , which is based solely on the X- 
ray brightness morphology. These are all cases in which 
Vikhlinin et al. (2002) classify a cluster as undisturbed, 
but more detailed X-ray (and in some cases radio) obser- 
vations show strong evidence of merging activity. Sec- 
ond, the clusters are classified according t o whether or 



not they have a central AGN, according to Mittal et al 
((2009i) . Additionally, we use the results of .Hudson et al 
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(|2010D . which divide the entire HIFLUGCS sample into 
strong, weak, and non-cool core clusters (SCC, WCC, 
and NCC) based on central cooling time, providing a 
third perspective from which to view the results. See 
Table |3] for each cluster's classification and the associ- 
ated references. 

4.2. Temperature Distributions 

Continuous temperature distributions for all clusters 
in our sample are shown in Figures H] - [71 The mea- 
sured temperature distribution is a combination of the 
intrinsic ICM temperature distribution and the ISF. To 
determine the relative contribution of the ISF to the 
shape of the temperature distributions, we directly com- 
pared the temperature distributions of the simulated 
clusters to clusters from the HIFLUGCS sample (Figure 
[8]). Least squares fitting to both the simulated cluster 
and HIFLUGCS distributions show they are both ap- 
proximately log-normal. Attempts to fit the distribu- 
tions with a single common function yielded similar 
values. Low temperature peaks are seen in many of the 
temperature distributions from the HIFLUGCS clusters, 
but we also see such secondary peaks in a few of the sim- 
ulations, indicating at least some of these peaks may be 
artifacts introduced in the MCMC, with the possible ex- 
ception of the binary clusters A399 and A401. However, 
the temperature distributions of the HIFLUGCS clusters 
tend to be slightly wider than that of similar simulated 
distributions. Thus, while the shape of the temperature 
distributions are clearly dominated by the ISF, we are 
capable of measuring extra broadening beyond the ISF. 
It is also worth noting that there are some variations 
in the temperature distributions from one cluster to the 
next; each cluster is unique at some level. 

4.3. Cluster Parameters 

We measured the median temperature, total emission 
measure, ^2500, and <TkT for each cluster in our sample 
(Table [3]). The median temperatures are also illustrated 
on the cluster temperature distributions (Figures H] - [7|). 

To evaluate the measurements of kTmed, we measure 
the median temperatures of the simulated isothermal 
clusters to test how well the known cluster tempera- 
ture can be recovered (Figure [9]). For clusters with 
kT > 5 keV, kTmed is biased slightly low. This is 
likely due in part to the high temperature side of these 
clusters' log-normal temperature distributions (see i i4.2p 
extending beyond the allowed temperature range for 
the smoothed particles, and the resulting asymmetry in 
the measured temperature distribution artificially low- 
ers kTmed- In addition, EPIC is not as sensitive to gas 
at such high temperatures, making it more difficult to 
constrain the temperature distributions. For real clus- 
ters, these two effects combine to lower the measured 
value of kTmed- Comparing kTm ed with the more tra- 
ditional te mperatures reported in fReiprich k, Bohringeii 
d^O^ and [Hudson etaLl (pOlOl ) (Figure©, we find rea- 
sonable agreement, although again for higher temper- 
atu re clusters kTmed ten d s to be slightly lower than 
the iReiprich fc Bohringeii (|2002D temperatures. In ad- 
di tion to the effect described abo ve, the temperatures 
of lReiprich fc B ohringeii ()2002[ ) and'HudsoneLalJ (|20100 
are obtained using completely different methods than 



this work, and may not represent the same physical tem- 
perature, and are also measured within larger apertures 
than r25oo, which is probably the main source of the dis- 
crepancy in the measured temperatures. There are sev- 
eral clusters (the seven nearest clusters and the Coma 
Cluster) for which ^2500 is completely beyond the spatial 
cuts described in iJ3 21 kTmed, cTkT, and EM for these 
clusters are therefore measured within a region smaller 
than r25oo, using emission from all available smoothed 
particles (within the spatial cuts). Excluding these clus- 
ters has little effect on the overall results (with the excep- 
tion of the slope of the EM — kTmed relation, see below), 
so they are included in the remainder of the analysis. 
Similarly, there are an additional 18 clusters for which 
^'2500 is only partially enclosed within the square region 
defined by the spatial cuts, either because r25oo is slightly 
larger than 10' or 15', or because the cluster observation 
was off center. However, in these cases some emission 
out to r25oo is included, and the amount of missing emis- 
sion area is at most ~ 5 — 10%. We therefore treat these 
clusters normally. 

From the perspective of our classification schemes, 
there are some notable differences in kTmed between 
typ es of clusters. A s found in other works, includ- 
ing i Chen et"al] ()2007[ ): iBurns et all (|2008l) : iMittal et al.l 
(2009), we find that relaxed, SCC, and AGN clusters 
(groups which largely overlap), all tend to have lower 
median temperatures, as well as smaller radii, than dis- 
turbed, NCC, and non-AGN clusters, as can be seen in 
Figures [inland El 

For a given cluster, akr quantifies any departure of 
the cluster ICM from isothermality. For all 62 clusters 
akT rules out isothermality by at least la, and in many 
cases much higher significance. Using the classification 
schemes introduced in §4.1) it can be seen from Figure [T2l 
that the akT appears to behave differently for different 
types of clusters. Relaxed, SCC, WCC, and AGN clus- 
ters are more likely to have narrower temperature dis- 
tributions (smaller akx), whereas disturbed, NCC, and 
non-AGN clusters are more likely to have wider tem- 
perature distributions. This effect is strongest for the 
disturbed/relaxed cluster classification. 

4.4. Scaling Relations 

We plot EM — kTmed, the analog of the Lx — T 
relation, in Figure 1131 Assuming a powerlaw relation 
logEM = alogkT -I- /3, a fit to all clusters u sing t he IDL 
MCMC fitting procedure linmix_err (Ke liyll2"007[ ) yields 
a = 2.69±0.24, /3 = 65.45±0.15, in agreement with other 
works. The slope is significantly steeper if only clusters 
with kTmed < 3 keV are included; a = 3.96 ± 0.98. As 
expected, relaxed, SCC, and AGN clusters tend to lie 
above and to the left of the best fit. If the eight largest 
clusters (for which r25oo is not visible) are excluded, the 
effect is to essentially remove most of the coolest clus- 
ters, which results in a shallower EM — kTmed slope of 
a = 2.18 ±0.29. 

Plots of akT ~ kTmed, cTkT " EM, akT - f2500 and are 
shown in Figures [141 [151 and [161 respectively. To inves- 
tigate the possibility of a correlation between akT and 
kTmed, the data was fit to two models, a linear relation 
and a constant. The best fit constant is akT ^1-38 keV 
with a of 333. The linear fit resulted in a best fit rela- 
tion of akT ~ 0.20kTmed+l-0S with a of 242. Though 
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neither the constant nor hnear model is a very good fit, 
the decrease in of Ax^ = 91 between the two fits hints 
at a correlation between akx and kTmed- The linear fit 
indicates the presence of intrinsic scatter of ~ 0.55 keV. 

5. DISCUSSION 

We have found that for relaxed, SCC, WCC, and AGN 
clusters, akx is more likely to be small, while for dis- 
turbed, NCC, and non-AGN clusters, (JkT is typically 
larger and may be correlated with kTmed- Apparently, 
the establishment of a strong cooling flow, often in con- 
junction with central AGN activity, leads to a reduction 
of both (TfcT and kTmed- As the gas cools, the higher 
temperature gas is removed, turned into lower temper- 
ature gas, down to s ome cutoff which scales as kTmed 
(jPeterson et alJl2003| ). likely due to AGN heating, which 
has the effect of lowering a^T- This would also have 
the effect of lowering kTmed- Nearly any disturbance to 
the ICM, generally from interaction with another clus- 
ter, usually disrupts the cooling flow, raising the median 
temperature. It also seems plausible that the merging 
of two clusters, with different temperature distributions, 
should result in a larger akT- Even for weaker interac- 
tions, where the clusters remain separate, the interaction 
may give rise to cold fronts, shocks, or at the very least 
a disruption to the cool cores, all of which could increase 
(TfeT- The presence of intrinsic scatter in UkT suggests the 
exact details of each cluster's merger, cooling, and AGN 
history probably has a signiflcant effect on auT- 

6. SUMMARY AND FUTURE WORK 

We have used the sm oothed particle inference method 
of UPeterson et al.l l200'^ to measure the temperature dis- 
tributions of the HIFLUGCS sample of galaxy clusters. 
The shape of the temperature distributions are found to 
be dominated by the isothermal spread function, which 
is approximately log-normal. The measured tempera- 
ture distributions tend to be somewhat broader than the 
ISF and have slight variations from cluster to cluster. 
We flnd the width of the distributions, CT/tr, is incon- 
sistent with isothermality for all 62 clusters. We flnd 
a slope for EM — kTmed relation that is steeper than 
predicted from self-similarity, but consistent with other 
works, a — 2.69±0.24, which increases to a = 3.96±0.98 
for low-temperature clusters. Separating the clusters ac- 
cording to various criteria, we flnd that the relaxed, SCC, 
and AGN groups largely overlap, as do the disturbed, 
NCC, and non-AGN cluster groups. The relaxed, SCC, 
WCC, and AGN clusters tend to have not only lower 
kTmed, as expected and found in other works, but also 
smaller akT than the disturbed, NCC, and non-AGN 
clusters. The akT — kTmed is described better with a 
linear flt, akT ~ ^-'^QkTmed + 1.08, than a constant, indi- 
cating that there may be a correlation between akT and 
kTmed- However, the intrinsic scatter in the akT mea- 
surements, ~ 0.55 keV, indicates that the exact history 
of each cluster may have a signiflcant impact on akT and 
there is still much to be learned. 

The upper end of the cluster temperature distribu- 
tions, particularly for massive clusters, is not as well con- 
strained as the lower temperature end. This is due partly 
to the fact that XMM-Newton does not have much sensi- 
tivity at energies above ^ IbkeV , and partly due to a lack 
of line emission at such high temperatures. However, it 



may be possible to somewhat improve temperature con- 
straints at higher temperatures by including hard X-ray 
data in the analysis, e.g. from Suzaku. The foremost is- 
sue to be addressed with future work is the large amount 
of scatter in the akT results, both to confirm or rule out 
a akT — kTmed correlation and to more accurately de- 
termine the amount of intrinsic scatter in akT - Towards 
this end, it would be useful to include more clusters in 
the sample, e.g. the 43 additional clusters in the HI- 
FLUGCS extended sample. The most useful avenue of 
research would be to compare our results with those from 
hydrodynamical simulations of cluster formation. Mea- 
surements of akT from such simulations would shed light 
on the amount of intrinsic scatter to be expected. They 
may also reveal how complicated cluster physics, such 
as conduction, AGN feedback, cooling flows, and cluster 
mergers, might effect akT- 
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ICM Temperature Distributions 



Table 1 

Observation Parameters 



Cluster 


Obs. ID 


Redshift 


■n-H 


PN 


Ml 


M2 


Rl 


R2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




uiiiiyuTUi 


0.00313 


1.75 


50.9 


58.6 


58.6 


55.3 


55.3 




U554d»U1U1 


0.00928 


4.91 


71.8 


103.8 


103.8 


108.7 


108.9 




Ulo215UlUl 


n r»i oori 

U. 012^9 


18.00 


19.0 


21.4 


21.9 


20.5 


23.4 


T?r~\ "D M A V 

r UixlN A A 


U4UUd2U1U1 


O.U04dU 


1.45 


68.7 


114.9 


114.9 


118.6 


118.6 




UUoU54UlUl 


n r»i /I ^ 

U.01d4d 


5.25 


26.4 


34.2 


34.2 




• • ■ 


A3581 


UoU4^5UoUl 


O.U2oOU 


4.30 


32.7 


53.3 


53.3 


„• ■ 

79.5 


79.5 


TV /TTjT'WTA 

MJ\ W 4 


UUyoUbUlUl 


U.02UU0 


1.86 


9.4 


14.1 


14.1 


14.5 


14.5 


A262 


UdU4/(5U1U1 


U.OluoU 


5.52 


8.3 


35.6 


35.7 


113.5 


113.5 


A400 


U4U4UiUiUi 


n no A Af\ 
U.U244U 


9.38 


20.5 


33.6 


34.4 


32.7 


32.7 


A1367 


r\r\r'-\ T/irn r\i 

UUd174U1U1 


0.U220U 


2.55 


17.6 


30.9 


31.5 






Zw54 


U5Uoz3U4Ul 


U.029U0 


16.68 


29.6 


45.3 


45.6 


48.8 


48.8 


TV /rT/"\S7'0 

MJ\ W o 


UoUUzlU/Ul 


n norrnn 
U.02 ([)() 


2.60 


16.5 


23. 1 


23. 1 


23.5 


23.5 


ASllOl 


U14/5UU1U1 


n nc^onn 


1.85 


80.0 


102.1 


102.8 


100.5 


100.5 


A3526 


U4Uu2UUlUl 


n m T /I n 


8.25 


84.3 


114.4 


114.4 


116.3 


116.3 


JAUooO 


U14/oUUzUl 


n nQ /I nn 
U.0o4y0 


24.00 


70.6 


79.9 


80.6 


106.1 


106.1 


A2052 


UlUy9zUlUl 


n no c: yi n 


2.90 


23.0 


30.2 


30.2 


25.8 


25.8 


A1060 


UzUuzoUiUi 


n n 1 ocn 
O.U12bU 


4.92 


28.9 


42.1 


42.1 


42.8 


42.8 


A1736 


UoUoziUzUi 


n n /I p; on 


5.36 


2.5 


10.5 


10.5 


12.0 


15.0 


A4038 


UzU44DUiUi 


n r\onr\r\ 
O.Uv5UOU 


1.55 


25.0 


29.4 


29.4 


29.7 


29.7 


MK Woo 


ni ririoom ni 

uiu9yduiui 


n n /I c nn 

0.04500 


3.15 


31.0 


32.0 


32.0 


31.4 


31.4 


A2634 


0002960101 


0.03139 


5.17 


6.0 


9.9 


9.9 






rL;X(JU42z 


U3U021U4U1 


n nonTn 

0.03970 


6.40 


31.0 


40.6 


40.6 


39.2 


39.3 


HYDA 


05042o0101 


n ncronn 

0.U5390 


4.86 


51.9 


82.4 


82.4 


95.6 


95.6 


AzUuo 


055U3d0101 


n no /I n ,1 

0.U3494 


2.92 


12.7 


23.8 


23.8 


26.5 


26.7 


A2597 


0147330101 


n no c: on 

0.U852U 


2.50 


48.3 


58.1 


58.2 


77.1 


77.1 


A2589 


0204180101 


0.04140 


4.39 


21.8 


32.8 


32.9 


32.0 


32.0 


A2147 


0505210d01 


n nocinn 

0.U350U 


3.29 


5.5 


10.7 


10.7 


11.6 


11.6 


A496 


U5UD2bU3Ul 


n noonn 

0.03290 


5.68 


39.5 


59.0 


59.1 


61.4 


61.4 


A O O'V/^ 

A3o7d 


01519U0101 


n n /I tr /^n 

0.04560 


5.01 


21.0 


21.0 


21.0 


27.2 


29.4 


A576 


U5U432U1U1 


n noonn 

0.03890 


5.69 


17.4 


33.8 


33.8 


16.8 


32.1 


A4059 


0109950101 


n n /I '7cn 

0.04750 


1.10 


12.6 


12.6 


12.6 


15.3 


15.3 


A2199 


UUUoUoUiUl 


n non i cr 

O.UoUlo 


0.84 


12.1 


13.9 


13.9 


13.8 


14.2 


A 1644 


UUiU42U2Ui 


n n /I '7Qn 


5.33 


11.1 


13.9 


14.6 




* * ' 


AO 133 


Ui44v5iUiUi 


n nc;^;£in 
O.UoddU 


1.60 


16.0 


22.1 


22.1 


22.1 


24.1 


A3395 


r\ ,1 n/~ir\ 1 r\or\i 

040U010301 


n ncrn/^n 

0.05060 


8.49 


22.3 


29.2 


29.2 


29.8 


0.0 


A3562 


U1U52d1oU1 


n n^innn 

0.U490U 


3.91 


34.6 


38.8 


39.2 




• • ■ 


A3112 


UuUoUoUlUl 


n nvc; o c: 
O.U 1020 


2.53 


63.2 


87.0 


87.1 


99.9 


100.2 


A3558 


01072d0101 


r\ r\ A onn 

0.04800 


3.63 


35.2 


43.6 


43.6 


37.5 


40.0 


A0119 


UoU5211(JUl 


r\ r\ A A on 

0.U442U 


3.10 


6.8 


8.8 


11.0 


12.2 


12.2 


A 1795 


UUy /oZUlUl 


n nco /I '7 
O.Ud24/ 


1.20 


30.4 


35.9 


36.0 


36.5 


36.5 


A2657 


U4UzlyUoUl 


n n/inon 
0.U4U2U 


5.27 


2.0 


22.8 


23.7 


24.3 


24.3 


A0085 


0Uo5140101 


n n cr c: nc; 

0.05506 


3.58 


9.0 


12.4 


12.4 






A3 158 


r\or\noi r\or\i 

UoUU21U2Ul 


n ncrn'7n 

0.U597U 


1.06 


10.4 


21.3 


21.3 


19.9 


19.9 


A3391 


UoUo21U4Ul 


n n T /in 
O.Uol4U 


5.42 


16.0 


26.0 


26.0 






A1650 


0U93200101 


n noo o ^ 

0.08384 


1.54 


30.6 


38.1 


38.1 


34.6 


35.6 


A2065 


0202080201 


0.07260 


2.84 


14.6 


20.6 


20.6 


20.7 


20.7 


A399 


Uii22DUiUi 


n n^zi c 1 
U.O ( lol 


10.58 


5.6 


13.1 


13.3 


9.0 


11.0 




r\<~>r\c orz r\-\ r\i 

020d850101 


n n cr /^n 

0.05560 


4.59 


47.8 


56.9 


56.9 






A3571 


r\noc;ncr r\or\i 

0U8d950201 


n non i n 

0.03910 


3.93 


12.2 


24.1 


24.1 


24.8 


24.8 


A1651 


0203020101 


0.08495 


1.71 


7.0 


10.7 


11.0 


12.9 


12.9 


A2256 


0141380201 


0.05810 


4.02 


10.8 


12.5 


12.5 


11.4 


11.8 


ZwC11215 


0300211401 


0.07500 


1.64 


16.9 


26.7 


26.7 


26.4 


26.4 


A401 


0112260301 


0.07366 


10.19 


8.1 


13.0 


13.0 


12.9 


12.9 


A2255 


0112260801 


0.08060 


2.51 


4.5 


12.2 


12.8 


11.5 


11.5 


A754 


0136740101 


0.05420 


4.59 


11.5 


13.8 


14.2 


15.6 


15.6 


A2204 


0306490201 


0.15216 


5.94 


13.6 


15.3 


15.3 


17.3 


17.3 


A2029 


0551780401 


0.07728 


3.07 


17.9 


24.5 


24.5 


44.0 


44.0 


A2142 


0111870301 


0.09090 


4.05 


4.5 


17.9 


17.9 






COMA 


0124711401 


0.02310 


0.89 


14.2 


17.4 


17.4 


19.1 


19.1 


A3266 


0105260901 


0.05890 


1.48 


15.0 


23.6 


23.8 






A478 


0109880101 


0.08810 


29.00 


59.0 


59.0 


65.4 


95.0 


103.6 


A2163 


0112230601 


0.20300 


12.27 


5.1 


10.3 


10.5 


7.8 


7.9 



Note. — (1) Cluster Name; (2) XMM-Newton Observation ID; (3) Optical red- 
shift, from the N AS A/IPACExtragalactic Database (NED); (4) Galactic hydrogen 
column density, in 10^" cm"'^, from lDickev Hz Lockmani (i 19901) , with the exception o f 
A478, NGC1550, and 2A0334, for which the value comes from lHudson et al.f l(20Tol '): 
(5)-{9) Net detector exposure times, in ks, after flare screening. 



Frank et al. 



Table 2 

Model Parameters 



Parameter Value Global 



ICM 

'n^ (0.?,nH,galA-'iriH,gal) yes 

log kT (keV), kTavg < 4 keV (-1.0, 1.0) no 

log kT (keV), kTavg > 4 keV (-1.0, 1.3) no 

Z (Z0)b (0.0,2.0) yes 

zned yes 

In a (arcsec) (0.5,5.5) no 

X (arcmin)'* (-10.0, 10.0) no 

X (arcmin) EPIC, large clusters'*,'' (—15.0,15.0) no 

y (arcmin) RGS<* (-2.3, 2.3) no 

y (arcmin) EPIC* (-10.0, 10.0) no 

y (arcmin) EPIC, large clusters'*, ° (—15.0, 15.0) no 



Galactic X-ray Background 

kT (keV) 016 yes 

Z (Zq)^ 1.0 yes 



Cosmic X-ray Background 

T L47 yes 

nH,gal yes 



EPIC Instrumental Background* 

Particle T -0.205 yes 

Ei 158 eV yes 



RGS Instrumental Background^ 

Particle F —0.35 yes 

Ei 500 eV yes 



Note. — Summary of model parameters and their allowed ranges 
(all priors are uniform). Values in parentheses represent the minimum 
and maximum values of the free parameters. All other parameters 
are fixed at the given values. Global parameters are the same for all 
smoothed particles in a given iteration. 

The absorbing column densities, riu g^i for each cluster, listed in 
Tabled 

*' M etallicitics arc with respect to solar metallicities of 
lAnde rs & Grcvcsse (1989). 

'^ The redshift is fixed for each cluster at the optical NED value. 

** The X and y coordinates are with respect to the cluster center. 

" For clusters which are large compared to the EPIC FOV (e.g. Coma), 

the X and y prior ranges were extended to 1 5' f rom the center to 

coincide with the spatial cuts to the data (see ^3.21 1. 

* All EPIC i nstrum ental background parameters are from 

lAndersson et al.l pOOTh . 

^ RGS instrumental background parameters are from [Peterson et al.l 
ll200l) . 

Table 3 

Cluster Properties 



Cluster 


fcT^ed (keV) 


o-feT (keV) 


EM (10*^^ cm"^) 


^•2500 (kpc) 


AGN 


Relaxed or Disturbed 


Cool-Core Status 


NGC4636 


0.72±0.06 


1.71±0.27 


0.33±0.05 


163.9±21.4 


Y 


Relaxed*'^ 


sec 


NGC5044 


1.22±0.02 


0.94±0.24 


4.90±0.41 


219.7±9.5 


Y 


Relaxed'' ■ 


sec 


NGC1550 


1.42±0.01 


0.82±0.20 


12.44±1.20 


239.6±8.8 


Y 


Relaxed* 


sec 


FORNAX 


1.48±0.03 


1.40±0.28 


0.90±0.10 


246.1±12.4 


Y 


Relaxed^ 


sec 


NGC507 


1.66±0.02 


1.48±0.29 


7.77±0.59 


262.9±10.1 


Y 


Disturbed®'^ 


see 


A3581 


2.00±0.01 


1.06±0.15 


34.31±2.10 


289.3±8.4 


Y 


Relaxed^ 


see 


MKW4 


2.10±0.02 


0.99±0.32 


12.51±1.06 


297.6±10.6 


Y 


Relaxed* 


sec 


A262 


2.36±0.01 


1.01±0.16 


21.24±1.21 


318.9±8.7 


Y 


Relaxed^ ■ 


see 


A400 


2.48±0.02 


1.35±0.19 


23.52±1.69 


326.6±12.1 


Y 


Disturbed*" 


Nce 


A1367 


2.56±0.02 


0.96±0.21 


22.22±1.16 


336.1±11.8 


N 


Disturbed** 


NCC 


Zw54 


2.56±0.02 


1.30±0.16 


43.64±1.98 


331.9±9.8 


Y 


Relaxed*^ 


wee 


MKW8 


2.74±0.03 


1.46±0.18 


27.07±2.15 


345.3±14.1 


Y 


Disturbed^ 


NCC 


2A0335 


2.76±0.01 


0.89±0.16 


662.01±25.88 


345.5±8.5 


Y 


Relaxed*^ 


sec 


ASllOl 


2.77±0.02 


0.83±0.22 


581.18±29.58 


342.4±10.5 


Y 


Relaxed*-* '2 


sec 


A3526 


2.88±0.01 


0.78±0.15 


19.10±0.93 


357.4±8.7 


Y 


Disturbed*5'2 


sec 


A1060 


2.98±0.04 


1.00±0.22 


17.04±1.35 


364.1±16.2 


N 


Relaxed*^ 


wee 
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Cluster 


fcT^ed (keV) 


o-fcT (keV) 


EM {lO'^'^cm-'^) 


^•2500 (kpc) 


AGN 


Relaxed or Disturbed 


Cool-Core Status 


A2052 


2.98±0.03 


0.86±0.23 


115.64±8.85 


360.7±15.3 


Y 


Relaxed^' 


sec 


A1736 


3.06±0.05 


1.80±0.23 


94.30±12.81 


364.1±22.5 


N 


Disturbed^ 


NCC 


A4038 


3.19±0.02 


0.91±0.14 


81.44±2.71 


375.6±8.3 


Y 


Tl T lis 

Relaxed^** 


WCC 


MKW3S 


3.27±0.05 


1.17±0.16 


172.40±15.81 


378.6±13.0 


Y 


Relaxed^-'''' 


sec 


HYDA 


3.49±0.02 


1.85±0.35 


311.31±14.54 


390.7±12.9 


Y 


1 i9n 9 
Relaxed^" ' 


sec 


A2634 


3.49±0.06 


1.56±0.27 


26.15±3.16 


400.6±29.7 


Y 


Disturbed^^ 


WCC 


EXO0422 


3.55±0.02 


0.98±0.14 


94.88±4.20 


397.6±9.3 


Y 


T-» 1 1 91 

Relaxed-^ 


sec 


A2147 


3.62±0.04 


2.95±0.41 


72.29±7.80 


402.8±20.2 


Y 


Disturbed^^ 


NCC 


A2597 


3.66±0.03 


1.14±0.44 


370.47±20.52 


395.4±16.4 


Y 


T-i 1 i99 9 

Relaxed^^ ' ^ 


sec 


A2063 


3.70±0.02 


1.23±0.15 


95.81±3.83 


407.6±10.6 


Y 


Relaxed^"^ ' 


WCC 


A3376 


3.78±0.03 


2.55±0.35 


74.49±4.49 


410.7±16.2 


Y 


Disturbed^^'-^'^ 


NCC 


A2589 


3.80±0.01 


0.97±0.14 


94.39±3.53 


413.0±9.2 


N 


Relaxed^S 


WCC 


A496 


3.81±0.02 


1.96±0.35 


140.83±6.31 


414.9±13.2 


Y 


T-» 1 T 97 

Relaxed-^ ' 


sec 


A576 


3.86±0.02 


2.09±0.33 


82.62±4.21 


416.8±10.8 


Y 


Disturbed'^" 


WCC 


A4059 


3.92±0.02 


1.87±0.31 


143.36±7.53 


419.1±10.7 


Y 


Relaxed" ' 
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Not e. — Measured clus ter parameters and classifi cations for our cluster sample. Quoted uncertainties are la errors. AGN classification 
is from lMittal et al.l (120091 '). Cool-core status is from [Hudson et al.l (12010') . 
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Figure 1. Ionization fractions of Fe K and Fe L as a function of temperature. These ionization fractions are used to weight the RGS and 
EPIC MCMC results as a function of temperature before combining. RGS smoothed particles are weighted by Fe L (solid, Wrgs)- Wepic 
is a piecewise function, such that EPIC smoothed particles which are within the RGS field of view are weighted by Fe K (dashed) and 
EPIC smoothed particles which do not overlap with RGS are weighted with the Fe K + Fe L ionization fractions (dot-dash). 
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Figure 3. ISF measurements for temperature priors of 0.1 < kT < 10.0 keV (asterisks) and 0.1 < kT < 19.5 keV (diamonds) as a function 
of the cluster temperature and tlie resultin g be st fi ts ( dotted and dashed lines correspond to the 0.1 < kT < 10.0 keV and 0.1 < kT < 19.5 
keV priors, respectively), as given in Eqns. Illl andll2l 
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Figure 4. Temperature distributions of clusters with — 10 keV prior for the smoothed particle temperatures. The vertical line marks 
feTmed ^■iid the dashed lines represent the Icr confidence interval of kTj^^^. Clusters are arranged in increasing order of kT^f,ii. 




Figure 5. Continuation of FigureU 
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Figure 6. Temperature distributions of clusters with — 20 keV prior for the smoothed particle temperatures. The vertical line marks 
fcT^gd aiid the dashed lines represent the la confidence interval of kT^f.^. Clusters are arranged in increasing order of kT^^^i. 
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Figure 7. Continuation of Figure [G] 
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Figure 8. Temperature distributions of the simulated isothermal clusters (black), with temperature distributions from the HIFLUGCS 
sample (the same as temperature distributions shown in Figures |4]-[7ll also shown for comparison (color). Each simulated cluster is shown 
with all clusters in our sample having kT^^,^ within 0.5 keV of the simulated cluster. The HIFLUGCS temperature distributions are 
vertically scaled to match the simulations and horizontally shifted such that the median temperatures are lined up. 




Fig ure 9. Comparison s of fcT^^^j from our sample with the average cluster temperatures reported in lHeiprich &: Bohringe3 I^OOZ) (left) 
and lHudson et aTl 1 12010 ) (middle), and for the simulated clusters, fcT^^jj with the simulated cluster temperature (right). The lines represent 
complete agreement between the two quantities (slope=l). 
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Figure 10. Histograms of fer^etj for each cluster subsample. Left: relaxed clusters (blue) and disturbed clusters (orange). Center: clusters 
witii (purple) and witfiout (green) central AGN activity. Right: SCC (blue), WCC (green), and NCC (red) clusters. 
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Figure 11. Histograms of r'2500 for each cluster subsample. Colors are the same as Figure [TOl 




Figure 12. Histograms of a/^x for each cluster subsample. Colors are the same as Figure [TUl 
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Figure 13. EM — kT scaling relation for all clusters in the sample, for each classification scheme. Stars represent the eight largest 
clusters for which r25oo is beyond the field of view. The powerlaw best fit to all clusters, logEM = alogkT^^^ + /3, is o = 2.69 ± 0.24, 
/3 = 65.45 it 0.15 (solid black). Including only clusters with kT^ed < SfceV, the best fit slope and intercept are a = 3.96 it 0.98 and 
/3 = 65.03 it 0.35, respectively (da shed). Excluding the eight largest clusters, the best fit is shallower, a = 2.18 it 0.29, /3 = 65.82 it 0.19. 
Colors are the same as Figure [lOl 
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Figure 14. kT^^d — <^kT scaling relation, shown with each classification scheme and the associated linear and constant fits, cf^t ~ 
0.20fcTmed + 1-08 (solid) and qj^ ~ 1.38 (dot-dash). Stars represent the eight largest clusters for which r2500 is beyond the field of view. 
Colors are the same as Figure [TOl 
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Figure 15. EM — cj-y scaling relation, shown with eac h cla ssification scheme. Stars represent the eight largest clusters for which r2500 
is beyond the field of view. Colors are the same as Figure 1101 
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Figure 16. r25oo — <^kT scaling relation, shown with eac h cl assification scheme. Stars represent the eight largest clusters for which r25oo 
is beyond the field of view. Colors are the same as Figure [TOl 



